arXiv:1507.00584vl [quant-ph] 2Jul2015 


Relation between the usual and the 
entanglement temperature, in a simple 
quantum system 


Alejandro Romanelli, Raul Donangelo, Andres Vallejo, 

Instituto de Fisica, Facultad de Ingeniena 
Universidad de la Republica 
C.C. 30, C.P. 11000, Montevideo, Uruguay 


Abstract 

We develop a thermodynamical theory to describe the behavior of the entanglement 
between a single two-level atom with a single mode of the electromagnetic held. The 
resonant Jaynes-Cummings model is used to study both the entanglement thermo¬ 
dynamics, in particular the entanglement temperature, and its connection with the 
average number of photons in the optical cavity. We hnd that this entanglement 
temperature has a strong dependence with the initial conditions of the atom. We 
show that the entanglement temperature between the photons and the atom dehned 
in this work is the same temperature obtained within the Jaynes-Cummings model 
at hnite temperature developed in the Thermo-Field Dynamics formalism. 
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1 Introduction 


Concepts such as thermodynamic equilibrium seem impossible to reconcile 
with the idea of isolated quantum systems since such systems follow unitary 
evolutions and do not reach a final stationary equilibrium state. Of course, 
a completely isolated quantum system is an idealization, constructed as a 
help to understand some phenomena displayed by real systems which may be 
regarded as approximately isolated. However, we recently [1,2,3,4] introduced 
the concept of temperature for an isolated quantum system which evolves in a 
composite Hilbert space. To do this we consider the quantum walk on the line 
(QW) (see [5] and references therein). The QW is a natural generalization of 
the classical random walk in the frame of quantum computation and quantum 
information processing and it receiving much attention recently [6,7,8,9]. 
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In our above mentioned works, we have developed a thermodynamic theory 
to describe the behavior of the entanglement between the coin and position 
degrees of freedom of the QW. Henceforth, we call “entanglement temper¬ 
ature” to the temperature associated to the entanglement between different 
degrees of freedom of an isolated quantum system. We have shown that, in 
spite of the evolution being unitary, in the QW a steady state is established 
after a Markovian transient stage. Those studies suggest that, if a quantum 
dynamics develops in a composite Hilbert space {i.e. the tensor product of sev¬ 
eral sub-spaces), then the behavior of an operator that belongs only to one of 
the sub-spaces may camouflage the unitary character of the global evolution. 
However it is not clear what is the relation between the usual temperature 
and the entanglement temperature. This question cannot be answered using 
the QW because it is an abstract mathematical model; to do this we need a 
real physical model where both temperatures, the usual and the entanglement 
temperature, emerge naturally. In order to answer it, we have here chosen one 
of the simplest and most interesting quantum models, the one known as the 
Jaynes-Cummings model (JCM) [10,11], that studies the interaction between 
radiation and matter. 

The JCM considers the interaction between a single two-level atom with a 
single mode of the electromagnetic held. The coupling between the atom and 
the held is characterized by a Rabi frequency, and a loss of excitation in the 
atom appears as a gain in excitation of the held. The collapse and the eventual 
revival of the Rabi oscillation, described by the analytical solution of the JCM, 
is a direct evidence of the quantum nature of radiation. The use of the JCM 
has permitted to elucidate basic properties of quantum entanglement as well 
as some aspects of the relationship between classical and quantum physics. 
Since it was proposed, the phenomenon has been of permanent interest in 
the quantum theory of interactions. About 30 years ago it was found that the 
model exhibits highly non-classic behavior, and the possibility of experimental 
realizations appeared. The relative simplicity of the JCM and its extensions 
has drawn much attention in the physics community and, more recently, in 
the field of the quantum computing [15,16]. 

Also in the 80’s, the Thermo Field Dynamics (TFD) formalism [12,13] was ap¬ 
plied to the JCM. The TFD is a method, developed in the 70’s by Takahashi 
and Umezawa [14], for describing Quantum Mechanical systems at hnite tem¬ 
perature. Using this method, it is possible to describe the statistical average of 
an observable at hnite temperature as a pure state expectation value. Thus, 
within the TFD formalism, one does not need to deal with a mixed state, 
which is a statistical ensemble of pure states at hnite temperature. In return 
for the above advantage, the TFD introduces the so-called tilde particles cor¬ 
responding to ordinary particles, thus doubling the dimension of the Hilbert 
space associated to the system. In the TFD method the ordinary particles and 
the introduced tilde particles represent the dynamical degrees of freedom and 
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the thermal degrees of freedom, respectively. 

In the present work we connect, within the JCM, the TFD thermodynam¬ 
ics with the entanglement thermodynamics presented in our previous works 
[1,2,3]. The paper is organized as follows. In the next section we review the 
usual JCM and study the photon thermodynamics in that model. In third 
section we develop the entanglement thermodynamics for the JCM and study 
its connection to the TFD-dehned temperature. Finally, in the last section we 
draw some conclusions. 


2 Jaynes-Cummings model 


We consider the ordinary JCM [11], composed by a single two-state atom in 
an optical cavity, interacting with a single quantized mode, with frequency u. 
The Hilbert space of the JCM has the form of a tensor product 

= ( 1 ) 


where the photon space, "Hn, is spanned by the unitary ortho normal vectors 
of the photon number state {|n)}, and the atom space, "Ha, is spanned by 
the two orthonormal quantum states {|e), |/)} that represent the excited and 
fundamental states of the atom, respectively. Note that the set {|n, e), \n, /)}, 
where \n,e) = \n)\e) and |n,/) = |n)|/), is an orthonormal base in the JCM 
Hilbert space. 


In this model, if the atom excitation frequency Ua is close to cj, then the system 
is near the resonance and it is possible to use the rotating wave approximation. 
In this case, and removing the field vacuum energy, the system Hamiltonian 
is 


H 


fc t ^ 

f\oj a'a + —oOaO'z 



( 2 ) 


where and a are the photon creation and annihilation operators respec¬ 
tively, and act on the photon number state \n). The radiation-matter coupling 
constant g is hxed by physical considerations such as the cavity volume and 
the atomic dipole moment. The raising and lowering operators are dehned by 


^- = \f){e\, (3) 

and the Pauli operator by 

= |e)(e| - |/)(/| = [a+,a_], (4) 
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and act on the atom states. Then the Hamiltonian, Eq. (2), is such that each 
photon creation is accompanied by an atomic de-excitation, and each photon 
annihilation by an atomic excitation. For a given photon number value n, Eq. 
( 2 ) has the eigenvalues 

E±(n) = ha;(n + ^) ± (5) 

where, for a specihc detuning parameter S = u — Ua, 

QriiS) = 1/52 {n + l)g‘^, (6) 

is the Rabi frequency. The corresponding eigenvectors, called “dressed states”, 
are given by 


n+} = cos 0 „|n, e) + sin9n\n + 1 , /), 

n—) = cos 6 '„|n-|-1,/) — sin 6 *„|n, e), (7) 


where 

tan(26'„) = ^ , ( 8 ) 

0 

Let us call |T(t)) the wave function of the JCM system. Its dynamics is given 
by the time-dependent Schrodinger equation and we will procure general so¬ 
lutions of the form 

00 

= E 

n=0 

where the coefficients and a~ are fixed by the initial conditions. In this work 
we consider the case of a separable atom-photon initial state. We assume that 
the initial state of the atom is arbitrary and that the initial photon number 
follows some initial coherent distribution. However, we must be careful to 
eliminate the | 0 ) |/) state in the initial condition because this state can not 
be built from the dressed states, see Eq. (7). More specihcally, we take initial 
conditions of the form 

|T( 0 )) =Aff^Cn\n) (^|e)cos| + |/)e*‘^sin|j -WColO) |/)e*^sin|, ( 10 ) 



iEj^ {n)t 


n+) + a^e 


iE_ (n)t 


|n—) 


(9) 


where the two parameters 7 G [ 0 , 7 r] and (p G [0,27r] dehne the initial state of 
the atom characterized by a point on the generalized Bloch’s sphere and M is 
a normalization constant. 


M 


1 

\/l- ICoPsmq 


( 11 ) 


Here we take {jCnP} as the initial coherent photon distribution inside of the 
cavity without atoms. Note that = 1, then n = is 
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the initial average number of photons. Using Eqs.(9),(10) and projecting on 
the base {|n, e), |n, /)}, it is straightforward to obtain the connection between 
{a+, a~} and the initial state of the atom and photon number, 


COS 9n COS ^ + Cn+1 siu 9, 


sin - I , 


a„ =J\f [ Cn+i cos sin — Cn sin cos ^ 


( 12 ) 


For this system, as is the case with most closed quantum systems, the prob¬ 
ability distribution does not converge in time. However we can use a natural 
notion of convergence in the quantum case, if we dehne the “limiting distri¬ 
bution” as the asymptotic limit of the average of the probability distributions 
at time 

^ ^ hs, 7 i*-PW*’ 

This definition captures the amount of time the system spends in each average 
state, and moreover, it corresponds to the natural concept of sampling from the 
system, since if one measures the system state at a random time chosen from 
the interval [0, t], the resulting distribution is exactly the average distribution. 
Applying the limit of the average to the previous dehnitions, we calculate the 
probability to obtain n photons independently of the state of the atom, i.e. 
P{n) = hm 4 ,_^.oo \ Jo • After a straightforward calculation, we obtain 

the hnal photon distribution in thermal equilibrium with the atom 

P(0) = Ioq pcos^ 6*0 -I- Ioq psin^ 6 *o, (14) 

F(n) = 9n + |a“psin^ 6 *„-|- |aj(j_;^psin^ 6 *n-i + |a“_ipcos^ 6 *„_i, (15) 

for n > 1. This hnal photon distribution depends both of the initial state of 
the atom as well as of the initial photon distribution, {|C*nP}- We can easily 
calculate the probability of hnding the atom in its excited (ground) state, 
independently of the number of photons in the system, i.e. 

^ oo oo 

Pe=lim-/ ^ |(n,e| T(t))|^dt = ^ |a+pcos^ 6 *„-F |a"psin^ 6 *„, (16) 

t Jo 


Pf = 1 - P.. 


(17) 


Note that, using Eq.(12), 

( c '^ cos^ 9n cos^ ^ -I- ^CnCn+i sm29n sin 7 cos v? -|- sin^ 9n sin^ 


(18) 


= (cjjsiii^^nCos^^ - ic„C„+isin2«„sin7COS43 + 7+, cos^#„siti^^j , 

(19) 
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then, the above probabilities are dependent both of the initial distribution of 
photons and of the initial state of the atom . 

In the asymptotic equilibrium, the average number of photons, (n) = nP{n), 

also can be calculated. Using Eqs.(14), (15), (18), (19), it is straightforward 
to show that and the average number of photons in the equilibrium satishes 

(n) = P/ + — (1 — e“”)sin^ 2^) ~ 2^) ’ 

In the case n 3> 1, it is easy to show that the difference An = (n) — n is 

An^ Pf- sin^ (21) 

so that, as it could be expected, the average number of photons does not 
change appreciably because of the presence of the atom. The probability of 
Ending the atom in its excited (ground) state in the resonant case is 

Pe = P, = 1 (22) 

independent of the atomic initial conditions. In this case <5 = 0 and 6n = 7^14: 
for all n, see Eqs.(8), (16) and (17). 

In the case n ~ 0 we have the vacuum of the optical cavity and Eq.(20) is 
reduced to 

(n) “ P,. (23) 

In this case, we can study the thermal effects in the JCM using the already 
mentioned TED formalism, according to which there is a correspondence be¬ 
tween the ensemble averages of statistical mechanics and the vacuum expec¬ 
tation values of quantum field theory [14]. This formalism was applied to the 
JCM in Refs. [12,13]. Ref. [13] obtains an explicit expression for the atomic 
state population as a function of temperature, namely 


Pe = 

g/3Sa;/2 

(24) 


Pf = 

g-/3&j/2 

(25) 

g/3/k<;/2 g—/3fitj/2 ' 


where [3 = {ksT) with ks the Boltzmann constant. 


3 Entanglement and temperature 


The unitary evolution of the JCM generates entanglement between the number 
of photons and the state of the atom. To characterize this entanglement we 
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start with the von Neumann entropy which is the quantum analogue of the 
Gibbs entropy 

5'Ar(p) =-tr(plogp). (26) 

where pit) = is the density matrix of the quantum system. Due to 

the unitary dynamics of the JCM the system remains in a pure state and this 
entropy vanishes. However, for these pure states the entanglement between 
the state of the atom and the number of photons can be quantihed by the 
associated von Neumann entropy for the reduced density operator that dehnes 
the entropy of entanglement 


S(p) = -tr(pclogpJ, 


(27) 


where here we dehne pc as 


Pc = lim - 

ti—>-00 f , 


pf 00 

•^0 n =0 


(28) 


and the partial trace is taken over the atom states. Using the wave function 
Eq. (9), together with Eq. (28), we obtain the reduced density operator 


pc 


^ Pe 0 \ 

V 0 ’ 


(29) 


Equation (29) is diagonal and we dehne the eigenvalues 


> 

+ 

III 

(30) 

III 

1 

< 

(31) 


For a given value of 5 7 ^ 0, these eigenvalues depend on the initial condi¬ 
tion of the atom through the angles 7 and p and also on the initial photon 
distribution, see Eqs.(16) and (17). 

In the resonant case, 5 = 0, the system stays in the degenerate case A_|_ = 
A_ = 1 / 2 , and in this case the entanglement entropy is maximum, S{p) = 1. 
If we want a more complete description of this equilibrium in the asymptotic 
limit, we need to know whether it is possible to associate a temperature to this 
entanglement entropy. In order to answer this question, we follow references 
[1,2,3]. According to the procedure described in them, it is necessary to connect 
the eigenvalues of pc with an unknown associated Hamiltonian operator He- To 
obtain this connection we use the quantum Brownian motion model of Ref. [18]. 
There one considers the system associated with the atomic degrees of freedom, 
(atomic states characterized by the density matrix pc) in thermal contact 
(entanglement) with the bath system associated to the photons. The number 
of photons can be thought as an external degree of freedom, in interaction 
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with the internal degree of freedom, the atom state. This picture is equivalent 
to the thermal contact between one system and its heat bath. In equilibrium 


Pq 0 , 


(32) 


should be satished. As a consequence, in the asymptotic regime the density 
operator is an explicit function of a time-independent Hamiltonian operator. 
If we note by {|*h±)} the set of eigenfunctions of the density matrix, the 
operators He and pc are both diagonal in this basis. Therefore the eigenvalues 
A± depend on the corresponding eigenvalues of He- We denote this set of 
eigenvalues by {e±}; they can be interpreted as the possible values of the 
entanglement energy. This interpretation agrees with the fact that A± is the 
probability that the system is in the eigenstate |*h±). The precise dependence 
between Aj. and e± = ±e is determined by the type of ensemble. In our case, 
this equilibrium corresponds to a quantum canonical ensemble. Therefore we 
propose that 


e±/3e 

= e/5^ + e-/5^’ 

which dehnes the entanglement temperature T = 1 /(/cb/ 9). Of course in Eq.(33) 
only the ratio e/T is well defined, however we choose to introduce this temper¬ 
ature as this concept strengthens the idea of asymptotic equilibrium between 
the photon distribution and the state of the atom. Note that while tempera¬ 
ture makes sense only in the mentioned equilibrium state, the entropy concept 
can be introduced without such a restriction. 


The probability that a state chosen at random from the ensemble {$+, *h_}, 
possesses an energy e is determined by the Boltzmann factor e~^^. Let us call 
'p'c the diagonal expression of the density operator pc, in the case treated in 
this paper p^ = pc, then 




1 

^ 0 \ 

pc = 

O 

> 

1 


y 0 e-^^j 


This operator is formally the same density operator that corresponds to an 
electron which possesses an intrinsic spin and a magnetic moment in an exter¬ 
nal magnetic held [17]. In general the Hilbert space of a quantum mechanical 
model factors as a tensor product Hsys ® Henv of the spaces describing the 
degrees of freedom of the system and the environment. The evolution of the 
system is determined by the reduced density operator that results from taking 
the trace over Henv to obtain Qgys = tienv{.p)- The simple models studied in 
Refs. [19,20,21] show how the correlations of a quantum system with other 
systems may cause one of its observable to behave in a classical manner. In 
this sense the fact that the partial trace over the photon number leads to a 
system effectively in thermal equilibrium, agrees with those previous results. 
Starting from Eq.(34) it is possible to build the thermodynamics for the JCM 



entanglement, where the partition fnnction of the system is then given by 


Z = e>‘‘ + e-^‘ = 2cos\i{l3e), 


(35) 


and the entanglement temperature is determined by 


T 



(36) 


When the system is in the resonance, A± = 1/2, it has maximum entropy and 
the entanglement temperature is also maximum, T = oo. It is important to 
point out that from the partition function Eq.(35) it is possible to develop 
all the entanglement thermodynamics, in other words, build the Hemholtz 
free energy, the internal energy and the entropy. As it should be expected, the 
entropy obtained from the partition function agrees with the previous Shannon 
expression, Eq.(40). Using Eqs.(30), (31) and (33), it is straightforward to show 
that 


Pe = 

e/36 

(37) 

1 

+ 

Pf = 

e-P 

(38) 

e/36 e-/36 • 


It is interesting to make the connection between Eqs. (37) and (38) with 
the equivalent expressions obtained within the TED formalism, Eqs. (24) and 
(25). From a straightforward comparison of these equations, we can determine 



Fig. 1. (a) Change in the average number of photons, (b) dimensionless entropy, 
and (c) inverse of the dimensionless entanglement temperature, l/T, as a function 
of 5. (d) Dependence of (n) on l/T. In all these calculations we took g = 0.001 and 
n = 0. 
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Fig. 2. (Color online) Five isothermal curves as function of the initial conditions, 
the dimensionless angles 7 and (p. Each curve is characterized by its thickness 
(and color). The values of /3 oc l/T, from the thickest to the thinnest line are: 
0.8 (purple), 0.7 (blue), 0.5 (green), 0.3 (yellow) and 0.1 (red). The isotherm cor¬ 
responding to /3 = 0 (T = 00 ) is situated between the two red isotherms. Here 
g = 0.001, n = 100,(5 = 0.01. 

and we conclude that the thermal vacuum temperature as defined in TFD 
and the entanglement temperature are essentially the same quantity, at least 
in the case of the JCM. The entanglement entropy, Eq.(27), can be expressed 
through the eigenvalues of pc as 

S (p) = -A+ log2 A+ - A_ logs • (40) 

Our results, for n ~ 0, are presented graphically in Fig. 1, which shows the 



Fig. 3. (Color online) The isothermal curves of Fig. 2 are shown on the Bloch sphere. 
The rotation axis of the figure is in the direction (1,0,1). The temperature of the 
entanglement between the atom and the photons is determined by the atomic initial 
condition. 


10 






behavior of (n) as a function of the resonant parameter 6. Figure 1(a) shows 
that for any value of 6, (n) < 1/2. Figure 1(b) shows that the entanglement 
entropy has its maximum value for the resonant case he. S' = 1 when <5 = 0. 
The maximum value of the entropy (maximum disorder) is achieved when 
the dimensionless entanglement temperature is T = oo, see Fig. 1(c). Under 
these conditions the atom excitation probability behaves as a classical Markov 
process. The relation between (n) and the inverse of the dimensionless entan¬ 
glement temperature, shown in Fig. 1(d), illustrates that, for large values of 
the entanglement temperature, a proportionality exists between (n) and 1/T. 
On the other hand for small values of T, (n) saturates to a constant value 0. 
These results are in accordance with Eq.(20). For n 3> 1, hgures 2 and 3 show 
the isotherms for the entanglement temperatnre as a fnnction of the atomic 
initial condition, Eq.(lO). In Fig. 2 the initial position is dehned throngh the 
angles 7 and (p and in Fig. 3 it is dehned throngh the position on the Bloch 
sphere. 


4 Conclusions 


We studied the unitary evolution of the JCM in a composite Hilbert space. In 
particular, we considered the entanglement between a single two-level atom 
and a photon distribution characterized by a single frequency. 

The JCM probability distribution does not converge with time. However, the 
asymptotic limit of the time-average of the probability distributions does con¬ 
verge. Thus the system establishes a stationary entanglement between the 
state of the atom and the photon distribution that allows to develop a ther¬ 
modynamic theory. The asymptotic reduced density operator is used to intro¬ 
duce the entanglement thermodynamic functions in the canonical equilibrium. 
These thermodynamic functions characterize the asymptotic entanglement. 
The JCM can be put into correspondence with a quantum walk on the line, 
where the internal degree of freedom (the state of the atom) is entangled with 
the position degree of freedom (the photon number). 

The difference between the average number of photons, with and without 
photon-atom interaction, depends on the entanglement temperature, and we 
show a connection between the entanglement temperature and the electro¬ 
magnetic held intensity in the optical cavity. 

We remark that in the case studied in this paper the entanglement temper¬ 
ature between the photons and the atom coincides with the temperature ob¬ 
tained within the TFD thermodynamics for the JCM. This implies that the 
entanglement temperature is amenable to direct measurements. We acknowl- 
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